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Discrete autonomous dynamical systems in dimension 1 can exhibit chaotic behavior, 
whereas the corresponding continuous evolution equations rule it out, and cannot even pos- 
qq I sess a nontrivial periodic solution. Therefore the passage from discrete to continuous equa- 

tions (and conversely) is all but harmless. We address this issue and evidence some caveats 

on the paradigmatic Verhulst logistic equation, investigating in particular the status and 

D ■ 

influence of the actual size of the unit time step in discrete modelings, rooted in well-known 
-(— > ' 

numerical analysis. 
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7 ' I. THE DISCRETE-TIME LOGISTIC EVOLUTION 

The logistic map f a (x) = ax(l — x) giving the celebrated recursion relation on the interval [0, 1] 
^H ! x n+1 = ax n (l - x n ) = f a {x n ) XQ 6 [0, 1] ffl6]l,4] (1) 

> ■ 

is one of the simplest example of discrete autonomous evolution leading to chaos. This nonlinear 

equation was introduced by Verhulst (a Belgian mathematician) in 1838 [ ] to take into account 

that a, the Malthus coefficient characterizing the growth of the population 
O ■ 
-i— > ' 

X n+ i = aX n , 

has to decrease when X n increases, due to resources limitation. The simplest way was to replace the 
constant rate a by a linear dependence in X n , matching the rate a at vanishing population, namely 
a(l — X n /M); the parameter M is then interpreted as being the maximum acceptable population, 

•i-H , 

^ currently known as the "carrying capacity" of the environment. Equation ^ is recovered through 

h ; 

the change of variable x n = X n /M. A very rich variety of dynamic behaviors is generated by 
Eq. (^Q), whose temporal structure is governed by the values of the control parameter a. Since 
the seminal reference [■], several studies of the asymptotic dynamics of (0) have been published, 
among which some very pedagogical ones are [£] and [*]. Let us only recall the most significant 
properties. 
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For a given, such that 1 < a < a\ = 3, the fixed point x* = 1 — 1/a is stable, globally attractive, 
therefore x n — > x a as n — > oo, irrespectively of the initial condition x provided it belongs to its 
basin of attraction ]0, 1[. In a\ = 3, a cycle of period 2 appears through a pitchfork bifurcation. 
Also called period-doubling bifurcation since it is associated with the destabilization of a fixed point 
x* a into a 2-cycle (or the destabilization of a 2 n -cycle into a 2 n+1 -cycle when it involves /„ instead 
of f a ), this generic bifurcation is characterized by the relation f ai (x ai ) = d x f(ai,x ai ) = —1 and 
the generic condition d ax f(a\,x* a ) ^ (denoting here the a-dependence on the same footing for 
the sake of clarity) [ 5 ]. The 2-cycle emerging in a\ remains stable and globally attractive in ]0, 1[ for 
any a < a 2 = l + \/6. More generally, there exists an increasing sequence (a,k)k of bifurcation values 
such that for au < a < ak+i, the asymptotic regime is a cycle of period 2 , which destabilizes in ctk+i 
through a pitchfork bifurcation of f£ . This sequence converges to Oqo ~ 3.5699 according to the 
scaling law Ooq — a^ ~ 5~ k with a universal rate 5 ~ 4.6692 [&]. The discrete evolution Q is actually 
a generic example exhibiting this so-called period-doubling scenario towards chaos, i.e. a normal 
form to which any one-parameter family experiencing such a scenario is conjugated [ 7 ]. In a = Oqo, 
a chaotic behavior arises, reflecting for a > a^ in a positive Lyapounov exponent (sensitivity to 
initial conditions) and mixing property (dynamic decorrelation of phase space regions). Chaotic 
regions in the a-space then intermingle in a highly complicated fashion (but now understood f 7 -]) 
with non chaotic regions where stable odd cycles rule the asymptotic dynamics. 

The conclusion, now acknowledged but striking at the time of publication of Ref. [ 2 ] and anyhow 
remarkable, is that a large variety of chaotic behaviors can be generated by a one-dimensional 
discrete evolution, with a seemingly harmless nonlinearity (smooth and simply quadratic). It 
showed that nonlinearities are never harmless when supplemented with a folding dynamics, here 
coming from the bell shape of the evolution map. 

II. CONTINUOUS-TIME COUNTERPART: A TRIVIAL DYNAMICS 

As it is impossible to give an analytical solution of (|T|). i.e. x n as an explicit function of n 
and Xq, and because we are interested in the asymptotic solution n — > oo (which gives a vanishing 
relative duration to the unit step n — > n + 1) it is appealing to deal with the corresponding 
continuous problem [ 9 ], which is straightforwardly solvable. To derive a continuous counterpart of 
(fT|. one subtracts x n to both sides of equation ((JJ) and identifies x n+ \ — x n with the differential of 
a continuous function of time y(t), which leads: 

f = f a (y)-y = y[a(i-y)-i], (2) 



whose analytical solution is easily obtained : 

m (a ~ 1)2/0 ( n\ 

V[t) ay + [a(l-y )-l]e-^)f W 

This solution is obviously regular with respect to t > for any value of a > 1 and, not surprisingly, 
tends to x* a when t — > oo. In contrast with this plain behavior, qualitatively insensitive to the 
value of a > 1, any attempt to solve (J2J) by discretization with a time step h = 1 will lead to the 
logistic evolution Q with its full richness of solutions as a is varied. On the other hand one expects 
that, for h small enough, one should approach the true solution (|3"|). How is it possible ? We have 
therefore to quantify what means "small enough" . 

III. INTERPRETATION OF DISCRETIZATION SCHEMES ASSOCIATED WITH THE 

LOGISTIC EQUATION 

Let us thus recall the behavior of the discretization schemes associated with (j^J^. Our aim is 
evidently not to get more knowledge about this equation, nor to device an accurate numerical res- 
olution, but rather to understand in this tractable and well-understood situation what is currently 
done to solve real problems when no straightforward solution is available. 

For a given time step h, the discretization scheme writes 

y(t + h)= y(t) + h{ay(t) [1 - y(t]) - y(t)} (4) 

A remarkable feature of the logistic equation is the possibility to rewrite this scheme as 

Y(t + h) = AY(t)(l-Y(t)), (5) 

with 

Y(t) = Xy(t) where \ = ^-J*— (6) 

involving the effective control parameter 

A(a,h) = l + h(a-l) (7) 

provided yo £ [0, 1/A] (note that A < 1 if h < 1). Obviously, the same phenomenology as for evo- 
lution ((J) will be observed. For instance, the inequality A < 3, required to obtain the convergence 
of © to the nontrivial fixed point Y£ = 1 — 1/A, means 

h < h c {a) = ^—!- = ~^-r (8) 

a — 1 a — 1 



Extending the reasoning to the subsequent bifurcations, one would observe a period-doubling sce- 
nario when the discretization step h increases, namely at values (hk)k with A{a,hk) = ctfe, i.e. 

h = ^ (9) 

a — 1 

Chaos arises for h > /ioo(o) = (ooo — l)/( a — !)• The bifurcation diagram as a function of h, at 

fixed a, would then be similar to the standard bifurcation diagram in a-space, up to a rescaling of 

the attracting sets by a factor of \(a, h), and a translation and rescaling of the bifurcation values 

(<Xfc = 1 + (a — l)/ifc). In particular, it is interesting to note that the sequence {hk)k follows the 

same universal scaling law hoo — hk ~ 5~ k or more precisely: 

h — h 
1+1 l ► 5 when i -> oo with 5 « 4.4669 (10) 

For illustration let us consider the case a = 3.1 (Figures 1, 2 and 3). The critical value of h 
is h c = (ai)/(a — 1) = 2/2.1 ~ 0.9524. For h > h c , one sees a 2-cycle, namely oscillations of 
the solution between the two (stable) fixed points of /^[/^(Y)]. The onset of the chaos is for 
h = h 00 = {a 00 - l)/(o - 1) = 2.5699/2.1 = 1.22376. 

IV. DISCUSSION: AN INTERPLAY BETWEEN TWO CHARACTERISTIC TIMES 

This simple study illustrates that the passage from continuous to discrete nonlinear equation 
is not insignificant: destabilization of the continuous time evolution, leading to cycles and even 
a spurious chaotic behavior, follows from an improper choice of the step of the discretization fii] 
or conversely an actual chaotic behavior can be suppressed by replacing a discrete model by its 
limiting continuous counterpart. 

Nevertheless, the passage from equation Q to (JSJ) by a simple scaling is exact only in the case 
of the quadratic family. We shall now investigate what remains in more general situations. Let / 
be a map, generating a discrete dynamical system x n+ \ = f(x n ) and having a stable fixed point 
x* (i.e. f(x*) = x* and |/'(x*)| < 1). The naive continuous counterpart writes dy/dt = f(y) — y. 
Linear stability analysis shows that x* is still a (at least locally) stable fixed point of the continuous 
dynamics since the linear growth rate of perturbations is negative: f'(x*) — 1 < 0. 

We might then consider the discrete scheme z n+ \ = z n + h[f(z n ) — z n ] for various values of the 
time step h. It is straightforward to show that this discretization scheme destabilizes for h > h c 
where 

h c = 2/[l - f'{x*)] (11) 



Indeed, the linear stability of x* breaks down when the modulus |1 + h(f'(x*) — 1)| overwhelms 1, 
which occurs for 1 + h(f'(x*) — 1) = —1. This relation yields the above value of h c and shows that 
the discrete scheme exhibits a period-doubling (pitchfork) bifurcation in h = h c (the additional 
generic condition for this bifurcation stated in Section 1 being also fulfilled, as can be directly 
checked) . 

The additional feature observed when the map f a depends on a control parameter a and exhibits 
a period-doubling in oi is that h c (a) crosses h = 1 in a = a±: for a > ax, f' a {x* a ) < — 1 and x* a is 
instable with respect to the initial discrete dynamics (h = 1) but is still a stable fixed point of the 
continuous dynamics, showing the inadequacy of the limiting continuous model dy/dt = f a (y) — y 
to capture the behavior of the discrete one x n+ \ = f a {x n )- It is to note that h c (a) decreases if a 
increases: the more stable is the fixed point (i.e. the larger \f' a {x* a ) — 1| with /a(x*) — 1 < 0), the 
smaller is the time-step range of validity of the discretization scheme (in a sense, the less stable is 
the discretization scheme). 

The qualitative differences, explicitly described in the previous sections, between the continuous- 
time and discrete-time versions of the logistic equation (and above in a more general framework) 
are not really surprising: a general claim assesses that a continuous-time dynamics requires a phase 
space of dimension at least 3 to develop a chaotic behavior [11]. In dimension 1 or 2, continuous 
trajectories behave as boundaries each for each other (trajectories of an autonomous continuous 
dynamic system cannot cross each other), which obviously prevents from chaos (and even nontrivial 
periodic solutions in dimension 1). But whereas it is thus straightforward to foresee the loss of 
chaotic and even periodic behavior when turning to the limiting continuous dynamics, is it possible 
to understand on physical grounds the existence of a critical value h c for the discretization time 
step hi The explanation lies in the comparison of the intrinsic time scale (s) of the dynamics with 
the chosen "time unit" h. 

The characteristic time of a continuous evolution, still denoted dy/dt = f(y) — y to avoid 
proliferation of new notations, can be estimated as r ~ 1/[1 — f'(x*)}. Indeed, a mere linearization 
of (J2J around the fixed point x* leads to: 

± [x (t) - x*] = [f'(x*) - l](x - x*) (12) 

hence the value of r. Destabilization of the discretization scheme occurs when h > h c = 2r. The 
stepwise updating, after each time step h, of the evolution law is too rough to properly control 
the discrete evolution and force it to follow closely all the relevant variations of the continuous 



trajectory. This is reminiscent of the Nyquist theorem fl£] for a periodic continuous evolution: the 
observation time step should be smaller than half the smallest period (or characteristic time) to 
properly sampling the continuous trajectory. 

It is to note that r or equivalently the critical value h c = 2r of the time step are intrinsic features 
of the dynamics, in the sense that they are invariant through conjugacy: for any diffeomorphism 
(j), f(y) and cfr^ 1 o / o 4>{y) (providing an equivalent modeling of the discrete model associated with 
/) or /(y) and y + 4>~ l [f o 4>{y) — 4>(y)] (providing an equivalent modeling of the continuous model 
associated with f(y) — y) will have the same critical value h c and the same characteristic time r. 

Let us carry further the comparison between the continuous evolution and its discretization, in 
order to understand the emergence of oscillations for h > 2r. The general continuous equation 
dy/dt = f(y) — y operates a fine tuning of the evolution rate dy/dt that is obviously not achieved 
by updating f(y) — y at times t n = nh. We have shown here that, near a stable fixed point, 
the resulting discrepancies lead to a bifurcation in the asymptotic dynamics, when h overwhelms 
the characteristic time of the evolution. To take a familiar example of such oscillations arising 
from a mismatch between two characteristic times, let us consider an heating/cooling device, 
able to measure the difference between the instantaneous room temperature and a prescribed 
one, and to monitor the appropriate energy supply or extraction, to compensate the measured 
difference. If the time h necessary for the device to actually deliver the required energy is longer 
that the characteristic time of temperature variations in the environment, the device will not 
balance the external temperature variations but rather, its ill-phased response will superimpose 
and the room temperature will suffer large oscillations. More generally, any ill-tuned homeostatic 
device, responding with a large lag h, will produce oscillations, and the result of Section 3 is the 
mathematical translation of this ubiquitous phenomenon. 

V. CONCLUSION 

In conclusion, we have presented an example showing explicitly the link between the validity of 
the discretization scheme with the dynamical (in)stability of the associated map for a unit step- 
size. Convertely, it enlights the specificity of the discrete dynamics, that cannot in general be 
understood, even qualitatively, from the behaviour of its continuous counterpart. In two or more 
dimensions, an additional problem arises: : the recursion relation is no more unique fi^]. The 
caveats illustrated in this paper are all the more relevant. 
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FIG. 1: Discretization of the logistic equation ^ with a = 3.1, using a time step h < h c (here h — 0.94 
whereas h c = 2/(a — 1) = 20/21 « 0.95J, see text, Section 3. Bold line: exact (continuous-time) solution of 
10). Stair step \fA{nh). x* a = 1 — 1/a. 




50 100 

FIG. 2: Same as Fig. 1 but with h — 0.96 > h c 
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FIG. 3: Discretization of the logistic equation 0) with a — 3.1, using a time step h — 3/(a — 1) = 1.424 
corresponding to the fully chaotic case A — A, see text, Section 3, and f]. 



